Diversity of Individual Mobility Patterns 
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Uncovering human mobility patterns is of fundamental importance to the understanding of epi- 
demic spreading, urban transportation and other socioeconomic dynamics embodying spatiality and 
human travel. The observed scaling laws for aggregated data require a theoretical explanation of 
their underlying mechanism. According to the direct travel diaries of volunteers, we show the absence 
of scaling properties in the displacement distribution at the individual level, which unfortunately 
provides a complete contrast to most inferences and assumptions in the literature. The aggregated 
displacement distribution follows a power law with an exponential cutoff, which is analytically ex- 
plained by the mixture nature of human travel under Maxwell-Boltzmann statistics. Our analysis 
provides an alternative way to bridge diverse patterns at the individual level and scaling laws at the 
population level. 
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I. INTRODUCTION 

Positioning systems in mobile phones and vehicles and 
Wi-Fi devices in laptop computers and personal digital 
assistants have made quantitative analyses of human mo- 
bility patterns possible PHI]. These analyses have a sig- 
nificant potential to reveal novel statistical regularities of 
human behavior, refine our understanding of the socioe- 
conomic dynamics embodying spatiality and human mo- 
bility k, 5] , and eventually contribute to controlling dis- 
ease designing transportation systems [l(| , locating 
facilities [ll[ , providing location- based services [HI, EH , 
and so on. 

Aggregated data from bank notes llj], mobile phones 
[l[ and onboard GPS measurements [2|] showed that the 
displacement distribution of human mobility, for both 
long-range travel and daily movements, approximately 
follows a power law. The scaling laws in long-range 
travel may result from the hierarchical organization of 
transportation systems [Hj], while the scaling laws in 
daily movements have recently been explained by the ex- 
ploration and preferential return mechanism (16J. This 
model [16| suggested that the displacement distribution 
at the individual level is power-law, which could also be 
considered to be a straightforward inference from the re- 
sults of aggregated data. 

Thus far, we still lack solid results about human mo- 
bility patterns at the individual level. Inferring individ- 
ual features from the aggregated data is very risky be- 
cause the scaling law for the population could be a mix- 
ture of many individuals with different statistics [13] • In 
addition, the aforementioned data are not sufficient to 
draw conclusions at the individual level. First, data such 
as GPS records from taxis and the trajectories of bank 



notes consist of many individual movements, but these 
individuals are not easy to be distinguished from each 
other. Second, data such as GPS records from mobile 
phones and the trajectories of bank notes could not ac- 
curately capture purposeful travels with explicit origins 
and destinations. In fact, the displacement between two 
activations of a mobile phone may be just a tiny portion 
of a purposeful trip or a combination of several sequential 
trips, while the displacement between two registrations of 
a bank note could be the result of a number of sequential 
trips made by different people. 
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Instead of using proxy data, we analyze the travel di- 
aries of hundreds of volunteers. Though the data set 
is small, it contains personal profiles and explicit posi- 
tions of origins and destinations, allowing quantitative 
and authentic analyses at the individual level. In con- 
trast to the scaling laws in aggregated data, people show 
diverse mobility patterns, and few of them display the 
scaling property. In fact, the trajectories of students and 
employees are dominated by trips connecting homes with 
schools and workplaces, respectively, while trips are dis- 
tributed more homogeneously among different locations 
for others such as retirees, homemakers, unemployed peo- 
ple, and so on. The aggregated displacement distribution 
follows a power law with and exponential cutoff, which 
can be analytically explained by the mixed nature of hu- 
man travel under Maxwell-Boltzmann statistics. In ad- 
dition, this theory predicts that the displacements using 
a single means of transportation will follow an exponen- 
tial distribution, which is also supported by the empirical 
data on taxi trips and air flights. 
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II. RESULTS 



A. Individual mobility patterns 



Our analysis of human mobility is based on a data 
set of 230 volunteers' six-week travel diaries in Frauen- 
feld, Switzerland [l8j]. This data set contains the volun- 
teers' personal information, including age, job and sex, 
and 36761 trip records. Each record includes the geo- 
graphic positions (longitude and latitude) of a trip's ori- 
gin and destination. By calculating the spherical distance 
between two trip endpoints from their longitudes and lat- 
itudes, we can obtain the length of each trip. 

We first measure the individual displacement distribu- 
tions from the data set. Figs. l(a)-l(c) show three typ- 
ical individuals' displacement distributions (Table SI 
presents all volunteers' displacement distributions), from 
which we cannot find any universal scaling properties. In- 
deed, when we use the Kolmogorov-Smirnov test [l9j to 
test whether the distributions fit the power law, we find 
that 87.8% of the individuals cannot pass the test (sta- 
tistical validation results are listed in Table S2). This 
result strongly suggests the absence of scaling laws in 
human travel at the individual level. 

To reveal the underlying structure of individual trips, 
we assign to each individual a mobility network, in which 
nodes denote locations visited by individuals, edges rep- 
resent the trips between nodes and edge weight is defined 
as the number of corresponding trips [20j. Figs. [Hd)-l(f) 
show three typical individuals' mobility networks (all net- 
works are presented in Table SI). As shown in Fig. Q] 
and Table SI, for most students and employees, their 
edge weights are highly heterogeneous. For each indi- 
vidual, we call the trip corresponding to the edge with 
the largest weight the dominate trip and define the dom- 
ination ratio d as the ratio of the weight of the dominate 
trip to the total weight. Fig. [5]reports the distribution of 
domination ratios for different groups of individuals, from 
which we can see that the student group has the largest d 
on average and that the employees' average domination 
ratio is smaller than that of the students but larger than 
that of the other group. 

The difference of d results from the fact that stu- 
dents and employees frequently travel between homes 
and schools/ workplaces in working days but retirees or 
homemakers do not have to do so. The peak values in 
the displacement distributions of students and employ- 
ees are thus usually determined by the lengths of their 
dominant trips. Because the lengths of dominant trips 
are not necessarily small, the displacement distribution 
for an individual is usually not right-skewed and is far 
different from a power law. In addition, the significant 
role of the dominant trip indicates that an individual's 
traveling process in general cannot be characterized by 
the Levy flight or truncated Levy flight. 
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FIG. 1. Individual mobility patterns, a-c. Displacement 
distributions for three typical individuals (a - a student, b 
- an employee, c - a retiree), where the peak values for the 
student and the employee result from the trips between two 
most frequently visited locations, d-f. Mobility networks for 
the three individuals, where the area of a node is proportional 
to its number of visits and the width of an edge is proportional 
to its weight. 
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FIG. 2. Distribution of the domination ratios, a. Population, 
b. Student group, c. Employee group, d. Others. S is the 
number of group members, and d is the average domination 
ratio. 
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FIG. 3. Displacement distribution P(r) of the aggregated 
data. The solid line indicates a power law with an exponential 
cutoff. 

B. Scaling property in aggregated data 

The aggregated displacement distribution of individu- 
als (see Fig. [3]) is well approximated by a power law with 
an exponential cutoff P(r) oc r~ 105 exp(— r/50), which is 
similar to those observed for bank notes [1J] and mobile 
phone users As shown above, this scaling property is 
not a simple combination of many analogous individuals. 
We assume that the total travel cost is C, the number of 
trips with cost Cj is Uj, and the same-cost trips are in- 
distinguishable, though they could be made by different 
individuals. Recalling Maxwell-Boltzmann statistics, the 
number of microstates is VL = N\ (EL "^) 1 where N is 
the total number of trips. According to the maximum 
entropy principle plj , the maximization of Q (under two 
constraints, ^ m — N and ^n^Ci = C) leads to the so- 
lution ni oc exp(— Ci/c), where c = C/N is the average 
travel cost. Denote the density of trips with cost c by 
P(c), then P(c) oc exp(— c/c). 

The travel cost is commonly approximated as the 
weighted sum c ks v/t + fj,m, where rj and fj, are two 
coefficients, and t and m are the costs involving time 
and money, respectively. Previous empirical studies have 
suggested that the monetary cost is approximately pro- 
portional to the travel distance as m w vr (22j, while 
the travel time approximately obeys a logarithmic form, 
t ?» cMnr + if) [2g, [13], where v, (f> and tp are coeffi- 
cients. The logarithmic relation results from the mix- 
ture of modes of transportation [25j . Apparently, peo- 
ple move faster when traveling longer distances: we walk 
from the classroom to the restaurant but take an air- 
plane from the US to China. Integrating the aforemen- 
tioned terms, we obtain the displacement distribution 
P(r) oc exp(— r/n), where /3 — rj(f>/c and n — c/ fiv. 

A direct corollary of Maxwell-Boltzmann statistics is 
that the displacement distribution should follow an ex- 
ponential form if it only accounts for trips from a sin- 



gle mode of transportation because in that case, c oc r. 
This corollary has found strongly supportive evidence 
from a number of empirical studies on disparate systems 
[2^ - [30j . Figure SI reports empirical distributions for 
taxi trajectories in Beijing [27], car trajectories in New 
York (downloaded from nhts.ornl.gov), bus trips in Shi- 
jiazhuang (collected by the authors) and air flights in the 
US (H)]. All distributions can be well characterized by 
exponential-like functions. 

III. DISCUSSION 

The general lessons that we learned from the present 
analysis could be used to refine our knowledge of hu- 
man mobility patterns. First, the displacement distri- 
butions for aggregated data usually display power-law 
decay with an exponential cutoff. Meanwhile, there are 
examples ranging from taxi trips to air flights in which 
the displacement distributions are exponential. In these 
examples, every displacement distribution is generated 
by trips involving a single mode of transportation, which 
corresponds to a linear relation between the travel cost 
and distance and eventually results in an exponential dis- 
placement distribution according to Maxwell-Boltzmann 
statistics. The present results suggest that the form 
(power law or exponential or other) of deterrence func- 
tion in the gravity law for human travel [3l| may be sensi- 
tive to the modes of transportation under consideration. 

This study warns researchers of the risk of inferring in- 
dividual behavioral patterns from aggregated statistics. 
Analogously, the temporal burstiness of human activities 
is widely observed, and the researchers are aware of the 
fact that the aggregated scaling laws could either be a 
combination of a number of individuals, each of whom 
displays scaling laws similar to the population (32j, or 
the result of a mixture of diverse individuals, most of 
whom exhibit far different statistical patterns than the 
population [33T - I351 ] . In comparison, such issues are less 
investigated for spatial burstiness. In particular, experi- 
mental analyses on individuals has rarely been reported. 
Determining whether the displacement distribution of an 
individual follows a power-law distribution will require 
further data and analysis. 

Many known mechanisms underlie the scaling laws 
of complex systems [3 6143 81 . including rich get richer 
[39l lit!, good get richer [4 ll lip , merging and regener- 
ation [431 ] . optimization [44l |45|, Hamiltonian dynamics 
[ifl , stability constraints |47| , and so on. The individual 
mobility model by Song et al. [l6[ is a typical example 
embodying the rich get richer mechanism. The maxi- 
mum entropy theory under Maxwell-Boltzmann statis- 
tics gives an unpretentious yet reasonable explanation for 
the emergent scaling from diverse individuals. In sum- 
mary, this work is complementary to known results on 
human mobility patterns and provides insights into how 
to bridge the gap between diverse individual statistics 
and aggregated regular patterns. 
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